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LONG WRITE-UP 



1 Introduction 



The harmonic polylogarithms (HPL) of Remiddi and Vermaseren introduced 
in [1] appear in many applications of high energy physics. They are a general- 
ization of the usual polylogarithms [2] and Nielsen polylogarithms [3] . Due to 
their definition through iterated integrations, they are useful for constructing 
solutions of differential equations, as they appear for example in the compu- 
tation of Feynman integrals. 

They are found in three-loop deeply inelastic splitting and coefficient functions 
[4,5,6,7], in two-loop massive vertex form factors [8,9,10,11,12,13,14,15], in 
two-loop Bhabha scattering [16,17,18,19,20,21,22], in multi-loop three-point 
and four-point functions [23,24,25,26,27,28,29,30,31,32,33], in 2- and 3-loop 
lepton g-2 [34,35], in the Higgs production and decay [36,37,38,39], heavy 
quark forward-backward asymmetry [40,41] and form factors [42], large- x limit 
of parton evolution [43] and various loop calculations [44,45,46] and in more 
formal developments [47]. 

HPLs appear in the expansion of hypergeometric functions, with real argu- 
ments for integer parameters [48,49,50,51] and imaginary arguments for hy- 
pergeometric functions with half-integer parameters [52,53]. HPLs of complex 
arguments also appear in two-loop massive master integrals [38,39] and in the 
anomaly contribution to the heavy quark form factors [13]. HPLs on the unit 
circle are related to the Clausen function and the generalized log-sine integrals 
[2,54,55,56,57,58]. 

The HPLs have been implemented in the computer algebra language FORM 
[59,60] and in the framework of GiNaC [61,62]. In a previous publication, we 
presented the Mathematica package HPL [63], which implemented HPLs for 
real arguments. The aim of this paper is to present its extension to complex 
arguments. 

The paper is articulated as follows. In the first part we shortly review the 
properties of the harmonic polylogarithms with the emphasis on the analytical 
properties in the complex plane. The second part describes the second version 
of the Mathematica implementation HPL of the harmonic polylogarithms that 
has been extended to treat arguments in the complex plane. 
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2 Review of the properties of the HPLs 



In this section, we review shortly the properties of the HPLs as described in 
[1,63]. 



2. 1 Definitions 



Harmonic polylogarithms (HPL) are defined [1] through recursive integration 
of so-called weight functions. The number of integrations is called the weight 
of the HPL. The usual weights are 



/i(*)= ' 



/o0r) = -, 

and the corresponding weight-one HPLs are 

X X 

if (1; x) = Jf 1 (t)dt = J — dt = - log(l - x), 



H(0;x) = log(x), 

X X 

H(-l;x)=ff. 1 (t)dt = J — dt = ]og(l + x). (2) 



HPLs of higher weights are then given by 



H( n Q-x) = — \og n x 
n\ 



H(a, ai,... tk ; x) = J f a {t)H(a h ...^ t) dt , (3) 



o 

where we use the notations 



n i = i,...,i, and a 1: ... :k = ^, a k - 

n 

A useful notation introduced in ref. [1] for harmonic polylogs with non-zero 
rightmost index is given by dropping the zeros in the vector a, and adding 1 to 



3 



the absolute value of the next non-zero index to the right for each dropped 0. 
This gives for example (3, —2) for (0, 0, 1, 0, —1). We can extend this notation 
to all index vectors by allowing zeros to take place in the rightmost position of 
the new index vector. This gives for example (3, —2, 0, 0) for (0, 0, 1, 0, —1, 0, 0). 
We will enclose index vectors in this notation in curly brackets and refer to 
it as the "m" -notation, as opposed to the "a" -notation. Some formulae or 
transformations are more easily expressed in the one or the other notation, 
therefore we keep both notations in parallel. 

For the extension to arguments in the complex plane it is convenient to define 
the linear combinations 



H(+;x) = H(l;x) + H(-l;x), 
H(-; x) = H(l; x) - H(-l; x), 
H(±, ai,..., fc ; x) = H(l, oi,... )fc ; t) ± H(-l, ai,... )fc ; t), (4) 

which correspond to introducing the weight functions 



We will refer to these two weights together with the weight f as the "±" 
weights as opposed to the "integer" weights /i, f and 

The formula for the derivative of the HPLs follows directly from the definition, 
-^F(a, ai,..., fc ; x) = f a (x)H(a li ... jk ; x). (6) 



2.2 Product identities and minimal set 



The product of two HPLs of same argument can be expressed as a sum of 
HPLs as follows, 

H(p 1 ,...,p Wl ;x)H(q 1 ,...,q W2 ;x) = H (p; x) H (q; x) = H ( r ; x )> ( 7 ) 

where p l±l q is the set of all arrangements of the elements of p and q such that 
the internal order of the elements of p and q is kept. For example one has for 
p = (a, b) and q = (y,z), 



H(a, b; x)H(y, z; x) = H(a, b, y, z; x) + H(a, y, b, z; x) 

+H(a, y, z, b; x) + H(y, a, b, z; x) 
+H(y,a,z,b;x) + H(y,z,a,b;x) . (8) 
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This identity is based on the general property 

t t t X\ t X2 

J dxi J dx 2 = J dx\ J dx 2 + J dx 2 J dx\ , (9) 



and is independent of the type of weights considered. 

As a consequence, the reduction to a minimal set for HPLs of ± weights is 
exactly equivalent to the reduction to the minimal set for integer weights. In 
our implementation of the ± weights, we use exactly the same procedure as 
for the integer weights [63], replacing the weight +1 by + and —1 by — . 



2.3 Analytical properties 



HPLs of integer weights have logarithmic divergences at 1,-1 and branch cuts 
on (— oo, —1), (1, oo) or both depending on whether they have 1 or —1 weights 
in their index vectors. HPLs with ± weights always have both branch cuts and 
have as many logarithmic singularities in both —1 and 1 as there are +s or 
— s before the first zero from the left of the weight vector. 

HPLs have logarithmic log n (0) divergences at and a branch cut on (— oo, 0) 
if they have n 0s at the right of the weight vector. 



2.4 Argument transformation 



This section presents the argument transformations for HPLs. We concentrate 
on the transformations for HPLs of ± weights. For integer weights we refer to 
refs. [1,63]. Since the argument of the HPLs are now complex, we have to treat 
the transformations in a more detailed way. We will take the argument to lie 
in the complex plane with the two cuts C_ = (— 00,— 1) and C + = (l,oo). 
We will further denote the upper half complex plane by C + and the lower 
complex plane by C_ . We consider the real axis on these cuts only as the limit 
of a complex argument approaching the real axis. So a "real argument with 
infinitesimal small imaginary part" will be addressed in this section exactly 
the same way as a argument in the upper complex plane. 

For HPLs on the real axis the argument transformations described in ref. [1] 
map the integration range to a range starting or ending at a point where a 
potential singularity lies. Therefore these transformations have to be applied 
to HPLs put into a form where the potentially divergent parts are explic- 
itly factored out in the form of HPLs of weight one, for which the analytic 
continuation is known. 
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Extending the HPLs to the complex plane allows to map the integration so that 
it avoids the divergent points for each HPLs, without having to put them in a 
factored form. In analogy to the "real" transformation which introduced HPLs 
evaluated at unity, the "complex" transformations introduce HPLs evaluated 
at i and —i. 



2.4-1 The transformation x — > — x 

We first consider the transformation x — ■» —x. It takes an argument from C + 
into C_ and conversely. Using the definition of the HPLs, 

r dx'f a (x')H(...,x') = - f dyf a (y)H(...,-y), (10) 
Jo Jo 

we see that the weight functions transform the following way 

fl ( X >) ^ -f^y) , / (x') -> f (y) , f^x') -> - fl (y) , 

UW) -> -/ + (y) , / (*') • • / (//)• (11) 

The transformation of HPLs without trailing 0s is easily done, we have for 
integer weights 

H(n 1 ,...,n w ]-x) = (-l) s H(-ni,...,-n w ;x) , (12) 

with s given by the number of non-zero weights (the length of the vector in 
the "m" notation) and for ± weights by 

H(s u s w ; -x) = (-l) s #(s!, s w ; x) , (13) 

with this time s given by the number of + weights only. So we see that HPL 
with ± weights and no trailing zeros are either odd or even, when the number 
of + weights is odd or even, respectively. 

The case of trailing zeros is more tricky. Consider x = re l<t> 

H (x) = log a; = log(r)+i0 , H (-x) = log(-rr) = log(r)+i(0±7r) , (14) 

with the + sign for argument in C_ and — for <C+. HPLs of higher weights 
with trailing zeros can be treated by extracting the trailing zeros using the 
product identities. 



6 



2.4-2 The transformation x — > 1/x 
We consider the transformation 

1 

x — — 

y 

for HPLs with ± weights. The results for integer weights can be obtained in 
the same way and are described in refs. [1,63]. The identities for weight one 
read 



H(0;x) = -H(0;y), 
H(+; x) = H(+; y) + m®(x), 

if (-; x) = H (-; y) + 2/7(0; y) + mB(x). (15) 

where Q(x) is +1 in the upper half complex plane and —1 in the lower half 
complex plane. For higher weight we proceed by induction. We take x to be 
in the lower complex plane, so that y is in the upper complex plane. We split 
the integration into two pieces, 



1/x 

H(s, s 2 ,...,kl V)= J d x'fs(x')H(s 2 ,..., k ; t) 






We have for the different cases s = and s = + and s — —, 



dt „ (\\ dt , , . , , 
- ¥ fo[ 1 )= T dt = f (t)dt, 

P- (?) = At {whji) = di (f^t - 1) = <->-« - 2 '°«> d <- 
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We only need to replace 

H ^ * ?) 

by its representation in terms of HPLs of argument t'. The case of x in the lower 
complex plane is treated exactly on the same way, separating the integration 
at x' = —i. 



2.4-3 The transformation x — > 

1 + x 

We consider the transformation 

X = 

1 + y 

for HPLs with ± weights. The results for the integer weights obtained on the 
same way and are described in [1,63]. The identities for weight one read 



H(0;x) = -H(+;y), 
H(+;x) = -H(0;y), 

if (-; x) = -if (0, y) + H (+, y) - ff (-; y) - 2 log 2. (18) 

For higher weight we proceed by induction. We split the integration into two 
pieces, 



1+x 



H(s,s 2 ... k ;y) = J dx'f s (x')H(s 2 ...k;t) 



l l + x \ 



J dx' + J dx' 



f s (x')H(s 2 _ k ;x') 



) 



H(s, s 2> ..., k ; i) -(j-j) J^iyfs (j^) H (s 2 ,..., k ; ^ 



(19) 



We have for the different cases s = 0, s = + and s — —, 
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2dt ' f t 1 -*')- 2dt ' - f (t > )At > 

j^f- (^) = « ^ = (/o(0 - U(t>) + /-(O) 

(20) 

Now we only need to replace 

l-t'\ 



H 8, 



,ki 



1 + f 



by its expansion in terms of HPLs of argument t' and perform the integration 
using the definition of the HPLs. 



2.4-4 The transformation x 2 — > x 

This transformation is not defined for ± weights. We mention it here since 
it will be used for values of HPLs at x = i in Section 2.5. The identities for 
weight one are 



ff(0; x 2 ) = log(x 2 ) = 2H(0;x), 

if (1; x 2 ) = - log(l - x 2 ) = ff (1; x) - if (-1; x) = ff (-; x). (21) 

The first identity holds on the complex plane except for the branch cut (— 00, 0) 
whereas the second equation holds on the entire complex plane except for 
(1, 00). For higher weights we use the relations, 



rx 2 fajJ 

H(0, m 2 ...,k\ x 2 ) = J —H(m 2 ,..., k ; x') 

= 2 ff*<- 

rx 2 (\ x > 

H{1, m 2 ..., k ] x 2 ) = - - — H(m 2 ,..., k ] x') 

= [" dt'f4t')H(m 2 t ;i' 2 ) (22) 

J 

recursively, where H(m 2 ^^t' 2 ) is expressed as HPLs of argument t', which 
are known in a recursive approach. 
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2.4-5 HPLs on the imaginary axis 



A nice property of the HPLs of weights + and — on the imaginary axis is that 
they are (provided they have no trailing zeros) either purely real or imaginary. 
Appending a + weight to an HPL that is real on the imaginary axis will change 
it into purely imaginary and vice versa. Appending a — weight doesn't change 
its real or complex nature. 

2. 5 Values at unity and i 
2.5.1 Values at unity 

HPLs of argument 1 are related to the Multiple Zeta Values (MZVs) described 
in the literature [64,65]. For positive weights, the value at unity is expressed 
in terms of usual MZVs and for general ms to colored MZVs. The relation can 
be found through induction and reads 

H({ mi _ k }- 1) = N( mi _ k )C(mi,..., k ), k > 1 (23) 
H({m};l) = ((m), m > (24) 

H({-m}; 1) = (1 -2 1 " m )C(m) m>0, (25) 

where the MZVs ( are defined by 

n A sgn(m>' , . 

11=1 12 = 1 «fe = l j=l % j 

The vector m is obtained from the vector m via 

m = (mi, sgn(mi)m 2 , . . . , sgn(mj_i)mi, . . . , sgn(m fe _i)m fe ). (27) 

The factor A^(mi v .. jfc ) is given by 

N(m lt ..., k ) = (-1)#(^<°). (28) 

The MZVs also form an algebra. Due to this fact, they can be expressed 
in terms of a few mathematical constants like powers of tt, ^-functions and 
polylogs at specified values. We list some of the identities of ref. [64,65] in 
Appendix B. For the implementation of the HPL at unity, we translated the 
tables of the FORM package harmpol .h [1] and their expansions for weight 7 
and 8 htable7.prc and htable8.prc to Mathematica. 

In these tables, there appear some constants that are not expressible through 
known constants like tt, ((n), log(2), or Li n (l/2). Using the different rela- 
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tions between the different MZVs, one can reduce the number of independent 
constants. The list of these independent constants are listed in Appendix B.2. 

For index vectors with leftmost Is, there appear divergences when the ar- 
gument approaches unity. Using the relations between HPLs and harmonic 
sums3, these divergences can be regulated in a proper way and expressed in 
terms of one divergent symbol 

oo -i 

5> = S(l,oo) = tf(l;l). 



Care has to be taken when taking the limit x — > 1 in expressions containing 
such divergences. Since we have, for example 



H(l,l;x) 



i) 2 



but 



(29) 



H(l, 1; 1) = 1 -H{1- If - Itt 2 . (30) 

we see that taking the limit x — ► 1 does not commute with the product identity. 
One should therefore express all the HPLs divergent in the limit x — > 1 in the 
same form (i.e. either with the logarithmic divergences factored out using the 
product identities, or with products of HPLs expanded in terms of a sum of 
HPLs of higher weight) and only then take the limit x — > 1. 



2.5.2 Values at i 

The argument transformations of Section 2.4 use the values of the HPLs with 
± weights at x — i. For the lower weights we have 



i)=log(2), (31) 

7i log(2)) , H(-, +,i)=i (-2C + ^tt log(2)) , (32) 

H(0,-,i) = ~, (33) 

2 

L., H {- 0,<) = ~-i7rlog(2). (34) 





111 


= Y ' 


H(+,-i) 


= i(2C 


H(0,+,i) 


= 2iC , 


H(+,0,i) 


= -2iC 



1 see section 5 of ref.[l]. 
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where C is Catalan's constant, 



c = T - 



To find the values at % for higher weights, one can use the following properties. 

• The product identities relate HPLs in % with HPLs of lower weights at %. 

• Some HPLs can be written in terms of more commonly known functions, 
whose values at x — i are known. 

• Using the argument transformation x 2 — > x one can relate the values of 
HPLs with only "-" and weights to the values of HPLs with integer weights 
at -1. 

• Since the "complex" transformations of section 2.4 involve HPLs at % and 
the "real" ones not, comparing them leads to relations between different 
HPLs at i. 

The number of relations is, however not sufficient to solve for all HPLs at %. 
A list of values of HPLs at i is given in Appendix A. 



3 The Mathematica implementation HPL 2 



In this section, we describe the extension of the Mathematica implementation 
HPL to complex arguments. Version 2 of HPL includes all features of the first 
version as well as the following new elements, 

• numerical evaluation of HPLs for complex arguments, 

• new weights + and — , 

• new function HPLArgTransf orm, 

• new integration procedure HPLInt. 

For completeness, we will not restrict the description to the new features of 
HPL but will also repeat the features that did not change. The package can be 
found at link [66] where installation instructions can be found. After installa- 
tion, the package can be loaded via 
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«HPL , 

*-*-*-*—*-* HPL 2.0 *-*-*-*-*-* 

Author: Daniel Maitre, University of Zurich 

Rules for minimal set loaded for weights: 2, 3, 4, 5, 6. 

Rules for minimal set for + - weights loaded for weights: 2, 3, 4, 5, 6. 

Table of MZVs loaded up to weight 6 

Table of values at I loaded up to weight 6 

$HPLFunctions gives a list of the functions of the package. 
$HPLOptions gives a list of the options of the package. 
More info in hep-ph/0507152 and at 
http : / /krone .physik . unizh . ch/~maitreda/HPL/ 



It should be loaded at the beginning of the Mathematica session. 



3. 1 New functions 



The package HPL defines the following new functions. 

• HPL[m,x] is the harmonic polylogarithm H(m;x), where m is a list repre- 
senting the index vector. For integer weights, we chose the "m" -notation as 
the standard notation. It is possible to give as the argument a vector in the 
"a" -notation, or even in a mix between the two notations, as no confusion 
is possible. Results will be displayed in the "m" -notation. It is also possible 
to mix ± and integer weights. 

HPL[{2, 1}, x] 
HPL[{1, 0, -1}, x] 
HPL [{plus, minus}, x] 
HPL [{plus, 3, minus}, x] 

HPL[{2, 1}, x] 

HPL[{1, -2}, x] 

HPL [{plus, minus}, x] 

HPL [{plus, 3, minus}, x] 



• HPLMtoA [m_List] and HPLAtoM [a_List] convert vectors from the "m"- to 
the "a"-, and from the "a"- to the "m"-notation respectively. Both can con- 
vert vectors which mix the two notations. 

HPLMtoA[{2, -3, 0}] 
HPLAtoM[{l, 0, 0, 1, 0, -1, 0}] 

{0, 1, 0, 0, -1, 0} 
{1, 3, -2, 0} 

HPLMtoA[{2, plus, minus, 0}] 
HPLAtoM[{0, 1, plus, 0, minus, 0}] 

{0, 1, plus, minus, 0} 

{2, plus, 0, minus, 0} 



• HPLpm21ml and HPL lml 2pm convert HPLs with ± weights into HPLs with 
integer weights and vice versa. 
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HPLpm21ml[HPL[{plus, minus}, t] ] 

-HPL[{-1, -1}, t] +HPL[{-1, 1}, t] -HPL[{1, -1}, t] +HPL[{1, 1}, t] 

HPLlml2pm[HPL[{l, -2}, x] ] 

1 1 
-— HPL [ {minus, 0, minus}, x] +— HPL [ {minus, 0, plus}, x] - 

1 1 

— HPL [ {plus, 0, minus}, x] + — HPL [{plus, 0, plus}, x] 



It is also possible to convert a function with integer weights into the corre- 
sponding sum of functions of ± weights, and vice versa. 



HPLpm21ml [ {plus , minus } 


, f] 




-f [-1, -1] +f [-1, 1] - 


f 


1, -1] +f [1, 1] 




HPLlml2pm[{l, -2}, f] 








1 

-—f [minus, 0, minus] 


+ 


1 1 

— f [minus, 0, plus] - — f [plus, 0, minus] 


1 

+ — f [plus, 0, plus] 



• HPLLogExtract extracts the singular behavior of HPLs using the product 
identities. The logarithmic divergences are extracted for 

integer weights in their argument at and 1, 
H — weights in their argument only at 0. 

The result is displayed as function of log(x), log(l — x) or H(l; x), H(0; x) 
depending on the option settings (see Section 3.3). 

HPLLogExtract [HPL [ { 1 , 2 } , x] ] 

HPL[{1}, x] HPL [{2}, x] -2 HPL [{2, 1}, x] 
HPLLogExtract [HPL [ {plus, minus, 0}, x] ] 

HPL[{0}, x] HPL[{plus, minus}, x] -HPL[{0, plus, minus}, x] -HPL[{plus, 0, minus}, x] 
$HPLAutoConvertToKnownFunctions = True ; 
HPLLogExtract [HPL [ { 1 , 2 } , x] ] 

-Log[l-x] PolyLog[2, x] - 

2 ( — Log[l - x] 2 Log[x] + Log[l - x] PolyLog[2, 1 - x] - PolyLog[3, 1 - x] + Zeta [3] J 

In the case of mixed integer and ± weights, only the divergences at are 
factored. 

• HPLConvertToKnownFunctions returns its argument with HPLs replaced by 
their representation in terms of more common functions, whenever possible. 



HPLLogExtract [HPL [{plus, minus, 0}, x] ] 




-HPL[{0, plus, minus}, x] - HPL [{plus, 0, 


minus}, x] + HPL [{plus, minus}, x] Log[x] 


HPLConvertToKnownFunctions [HPL [ { 2 , 2 } , 


*]] 


1 7 

— PolyLog[2, x] z -2PolyLog[2, 2, x] 




HPLConvertToKnownFunctions [HPL [ {plus} , 


x]] 







It is only needed if the option $HPLAutoConvertToKnownFunctions is set 
to False (see Section 3.3). 
• HPLProductExpand returns the value obtained by replacing products of 
HPLs of same argument by their representation as a linear combination 
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of HPLs, as presented in Section 2.2. 



HPLProductExpand[HPL[{2, 1}, t] HPL[{1, 0}, t] ] 

HPL[{1, 2, 2}, t] +2HPL[{1, 3, 1}, t] +2HPL[{2, 1, 2}, t] + 
HPL[{2, 2, 1}, t] +HPL[{1, 2, 1, 0}, t] +3HPL[{2, 1, 1, 0}, t] 

HPLProductExpand[HPL[{l, 0}, t] HPL[{-1}, t] HPL[{1, -1}, t] ] 

2HPL[{-1, 1, 1, -2}, t] +HPL[{-1, 1, 2, -1}, t] +HPL[{1, -2, 1, -1}, t] + 
3HPL[{1, -1, 1, -2}, t] +HPL[{1, -1, 2, -1}, t] +4HPL[{1, 1, -2, -1}, t] + 
4HPL[{1, 1, -1, -2}, t] +2HPL[{1, 2, -1, -1}, t] +HPL[{-1, 1, -1, 1, 0}, t] + 
2HPL[{-1, 1, 1, -1, 0}, t] +2HPL[{1, -1, -1, 1, 0}, t] + 
3HPL[{1, -1, 1, -1, 0}, t] +4HPL[{1, 1, -1, -1, 0}, t] 

HPLProductExpand[HPL[{plus, minus}, t] HPL[{1, -1}, t] ] 

HPL[{1, -1, plus, minus}, t] +HPL[{1, plus, -1, minus}, t] +HPL[{1, plus, minus, -1}, t] + 
HPL[{plus, 1, -1, minus}, t] +HPL[{plus, 1, minus, -1}, t] + HPL [ {plus, minus, 1, -1}, t] 



In order to expand all products, HPLProductExpand expands its argument 
(using Expand), so that terms of the form 



H(...;x)(H(...;x) + H(...;x) + ...). 



are also replaced. For large expressions this might be time consuming, one 
should in this case first collect the products of HPLs and apply HPLProductExpand 
only to the products and not to the whole expression. 



HPLProductExpand[HPL[{2}, t] (HPL[{1, 0}, t] +HPL[{-1}, t] ) ] 

HPL[{-2, 1}, t] +HPL[{-1, 2}, t] +2HPL[{1, 3}, t] + 
HPL[{2, -1}, t] +HPL[{2, 2}, t] +HPL[{1, 2, 0}, t] +2HPL[{2, 1, 0}, t] 



Since the product properties are general, the function works for any vectors 
of symbolic weights. 



HPLProductExpand[HPL[{a, b}, x] HPL[{y, w, z}, x] ] 

HPL [{a, b, y, w, z}, x] +HPL[{a, y, b, w, z}, x] + 
HPL [{a, y, w, b, z}, x] +HPL[{a, y, w, z, b}, x] + 

HPL[{y, a, b, w, z}, x] +HPL[{y, a, w, b, z}, x] +HPL[{y, a, w, z, b}, x] + 
HPL[{y, w, a, b, z}, x] +HPL[{y, w, a, z, b}, x] +HPL[{y, w, z, a, b}, x] 



• HPLConvertToSimplerArgument returns its argument with HPLs of related 
arguments replaced by their expansion as a sum of HPLs of simpler ar- 
guments. Which transformations are implemented depends on the type of 
weights. They are: 

integer weights —x, x 2 , 1 — x, 1/x, x/(x — 1) and (1 — x)/(l + x) 
± weights — x, 1/x, and (1 — x)/(l + x). 
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HPLConvertToSimplerArgument [HPL [ { 2 , 1 } , - x] ] 

HPL[{-2, -1}, x] 

HPLConvertToSimplerArgument [HPL [ { 2 , 1 } , 1-x]] 

-HPL[{1, 0, 0}, x] + Zeta[3] 

1-x 

HPLConvertToSimplerArgument [HPL[{1, 1}, 11 

1 + x 

-HPL[{-2}, x] +HPL[{-1, -1}, x] -HPL[{-1, 0}, x] + 

HPL[{0, 0}, x] -HPL[{-1}, x] Log [2] +HPL[{0}, x] Log [2] + L ° g j; 2] 
HPLConvertToSimplerArgument [HPL [{2}, (x) / (x - 1 ) ] ] 

-HPL[{2}, x] -HPL[{1, 1}, x] 

HPLConvertToSimplerArgument [HPL [ {plus, minus}, (1-x) / (1+x)]] 

7T 2 

— + HPL [ {0, 0} , x] + HPL [ {0, minus} , x] - HPL [ {0, plus}, x] + 2 HPL [ {0} , x] Log [2] 
6 



The transformations do not work for mixed integer and ± weights. 
HPLArgTransf orm [x , r ,5 ] returns its first argument expressed in terms 
of HPLs of the transformed argument. It does the same job as HPLConvertToSimplerArgument, 
but the transformation to be performed is not fixed by the form of the ar- 
gument of the HPL but is specified as an argument of HPLArgTransf orm. 



HPLArgTransform[HPL[{2, -1}, X], x -x] 
HPLArgTrans form [HPL [ {plus, 0, minus}, X], x-»-x] 

HPL[{-2, 1}, -X] 

-HPL [{plus, 0, minus}, -X] 

HPLArgTransform[HPL[{2, 1}, X], x -> Vx] 

2 (HPL[{-2, -1}, Vx] -HPL[{-2, 1}, Vx 7 ] -HPL[{2, -1}, Vx~] +HPL[{2, 1}, Vx] 
HPLArgTransform[HPL[{l, 1, 2}, X], x -> 1 - x] 

4 1 

^— + — 7T 2 HPL [{0, 0}, 1-X] + HPL [{3, 0}, 1-X] + 2 HPL [ {0}, 1-X] Zeta[3] 
30 6 



In addition one can specify whether the argument is assumed to lie in the 
upper (8—+1) or lower (5 = —1) complex plane, or whether to leave it 
purely symbolic. 
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HPLArgTransform[HPL[{2, 2}, X], x -» x / (x - 1) ] 

hpl[{2, 2}, +HPL[{1, 1, 2}, — x 



X 

hpl[{2, i, i}, _ 1 X +X ] +HPL[{1, 1, 1, 1}, t^t] 

HPLArgTransform[HPL[{2}, X] , x -» (1 - x) / (1 + x) , 5] 

1-X, „„ r , , „, 1-X n , , 1-X 



HPL 



[{-1, -1}, y^] +HPL[{-1, 0}, TTY ] -HPL[{1, -1}, 



1 + X 



HPL[{1, 0}, -i— |] + HPL[{-1}, 1— |] Log[2] + HPl[{1}, 1— |] Log[2] 
HPLArgTransform[HPL[{plus, 0, minus}, X] , x -» (1 - x) / (1 + x) , 5] 

- A ^-^- 2 HPL[ { } , l^f]-HPL[(0, plus, 0}, l^f]- 

1 — X 1 — X 

HPL[{0, plus, minus}, — 1 +HPL[{0, plus, plus}, 



1+x j l L ' * J ' 1 + X J 

2HPL[{0, plus}, Log[2] - 7Zet 4 a[3] + J-izeta[3, -1] - ^ iZeta[3, -|] 

HPLArgTransform[HPL[{2, -2}, X], x -» 1 / x, <5i ] 

+HPL[{-4}, 1] - i-^ HPL[{2}, 1] - i-^ HPL[{0, 0}, A-] + HPL[{2, -2}, ^] - 

HPL [{2, 0, 0}, ^] -HPL[{0, 0, 0, 0}, ] - -L i tt 3 HPL [ { } , ^ ] 6 1 + -| HPL [ { } , ^] Zeta[3] 
HPLArgTrans form [HPL [ {minus, 0, plus}, X], x-»l/x, <5 2 ] 

-2HPL[{0, 0, plus}, -HPL[{minus, 0, plus}, - ^ 1 ^ 6 2 " 

2 i ttHPL [ {0, 0} , — ] 5 2 - i tt HPL [ {minus, 0} , — ] 5 2 + Ji 2 HPL [ {0}, — ] 5\ + 

-i- tt 2 HPL[{minus}, ^ ] &\ + ±- i tt 3 6\ - i 6 2 Zeta[3, |] + ^i<5 2 Zeta[3, -|] 

The transformations available are the same as for HPLConvertToSimplerArgument, 
with the same restrictions on the weight vector. 

HPLArgTransform[HPL[{2, -1}, X], x-»Vx] 

HPLArgTransf orm: :WrongWeight : 
The argument transformation X -Jx attempted is not defined for the weight vector {2, -1} 

HPL[{2, -1}, X] 

HPLArgTransform[HPL[{l, -1, 2}, X], x -» 1 - x] 

HPLArgTransf orm: :WrongWeight : 
The argument transformation X^l-X attempted is not defined for the weight vector {1, -1, 2} 

HPL[{1, -1, 2}, X] 

HPLArgTransform[HPL[{2, -1}, X] , x -> (x) / (x - 1) ] 

HPLArgTransf orm: : WrongWeight : 

The argument transformation — - attempted is not defined for the weight vector {2, -1} 
HPL[{2, -1}, X] 

The rule r can be written using any symbol, which makes the code more 
intuitive. The symbol in the rule does not have to match the arguments of 
the HPLs. 

HPLArgTransform[HPL[{2, 3}, - Vx ] , t -> -t] 

HPL[{-2, -3}, Vx"] 

HPLArgTransform[HPL[{2, 3}, 9/10], x 1 - x] 

^/T 4 HPL[{1}, J_] -1^ H PL[{1, 2}, A-] - 

HPL[{1, 2, 1, 0}, -L] - -L/t 2 Zeta[3] +HPl[{1, 0}, -^} Zeta[3] + 9 Zet 2 a[5] 
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• HPLReduceToMinimalSet returns its argument with the HPLs projected to 
the minimal set, as described in Section 2.2. 

HPLReduceToMinimalSet [HPL [{1, 2 , 3}, x] ] 

HPL [ {2}, x] 3 - 4 HPL [ {3} , x] HPL [ {2, 1} , x] - HPL [ {2} , x] HPL [ {2, 2} , x] + 
HPL[{1}, x] HPL[{2, 3}, x] +2 HPL [{2, 1, 3}, x] -2HPL[{2, 3, 1}, x] 

HPLReduceToMinimalSet [HPL [ {minus, minus, 0}, x] ] 

1 7 

— (HPL[{0), x] HPL[{minus}, x] - 

2 (HPL [ {minus } , x] HPL [ { , minus } , x] - HPL [ { , minus, minus } , x] ) ) 

• HPLAnalyt i cCont inuat ion [x , HPLAnalyt i cCont inuat ionRegion->region] 

assumes that the arguments of the HPLs are real and returns its argu- 
ment x with HPLs replaced by their analytic continuation. The arguments 
of the HPLs are taken to belong to the interval specified by the option 
HPLAnalyticContinuationRegion, for which region can be either 
minftoml the interval — oo to —1 
mltoO the interval —1 to 
onetoinf the interval 1 to oo. 

The HPLs are replaced by their representation in terms of HPLs of argument 
in the interval (0, 1). The choice of the side of the branch cut from which 
the argument is approached is set by the option 

HPLAnalyt i cCont inuat ionS ign 
which can take values —1, 1 or any symbol. 

HPLAnalyticContinuation [HPL [ { } , y] , 
AnalyticContinuationRegion -» mltoO, AnalyticContinuationSign -» 1] 

±71+ HPL [ {0} , -y] 

Simplify[HPLAnalyticContinuation[HPL[{l, -1}, y] , 

AnalyticContinuationRegion -» onetoinf, AnalyticContinuationSign -» -1] ] 

_ZL_ +HPL[{-2}, — ] -HPL[{0, 0}, — ] +HPL[{1, -1}, — ] -HPL[{1, 0}, — ] -i7TLog[2] 

If the option HPLAnalyticContinuationRegion is omitted, and if the ar- 
gument is numerical, HPLAnalyticContinuation will automatically use the 
appropriate setting. If the option HPLAnalyt icCont inuat ionS ign is omit- 
ted, HPLAnalyticContinuation will use the value stored in the variable 
$HPLAnalyticContinuationSign which is set by default to 10- 



HPLAnalyt icContinuation[ HPL [ {1} , — ] , AnalyticContinuationSign 1] 

-i7T + HPL[{0}, y] +HPL[{1}, y] 

5 

Simplify [HPLAnalyticContinuation [HPL[{1, 0}, -— ]]] 

^-HPL[{-2}, -|] +i (tt + i HPL[{0}, (-HPL[{-1}, -|] +HPL[{0}, -| ] ) + HPL[{0, 0}, |-] 

It is to be noted that the Mathematica conventions for the analytic continu- 
ation are not always the same as that of the HPL package. This is illustrated 
by the following example 



2 This is the same convention as ref. [67], but opposite to that of ref. [62] 
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HPL[{1}, 1.9] 

0.105361 + 3.14159 i 
-Log[l - 1.9] 

0.105361 - 3.14159 i 

Here HPL takes the argument of the HPL to have an infinitesimal positive 
imaginary part, whereas Mathematica take the argument of the logarithm 
to have a positive imaginary part. Since the substitution of HPLs through 
more common functions has precedence over the analytic continuation, the 
option 

$HPLAutoConvertToKnownFunctions 

can interfere with the analytic continuation. 

$HPLAutoConvertToKnownFunctions = True; 
HPLAnalyticContinuation [ HPL [ { } , t ] , 
AnalyticContinuationSign -> <5, AnalyticContinuationRegion -> mltoO] 

Log[t] 

$HPLAutoConvertToKnownFunctions = False; 

HPLConvertToKnownFunctions [HPLAnalyticContinuation [HPL [ { } , t ] , 
AnalyticContinuationSign-* <5, AnalyticContinuationRegion-* mltoO] ] 

i 7T 6 + Log [ -t ] 

This example shows that with the option 
$HPLAutoConvertToKnownFunctions 

set to True we lose control over the sign of the imaginary part (as in the 
first case it will now depend on Mathematica's conventions). 
• HPLInt [x,t] returns a primitive of x with respect to t. The cases working 
are for k > 
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HPLInt [y + y / (1 - y) A 3 HPL [ {plus, minus}, y] , y] 

y 2 y (-5 + 6 y) (2 - 3 y) HPL [ {minus}, y] 
2 + 8(-l+y) 2 + 4 (-1 + y) 2 + 

5 HPTTfnluc.) vi i (~l + 2y + 3y 2 ) HPL [{plus, minus}, y] 
— HPL [ {plus} , y] + 8 (_l + y) 2 

HPLInt f HPL [{plus, minus}, y] , yl 

y(i + y) 2 

/ y HPL [{minus}, y] 1 nnT 
c ( TTy + — y HPL[{plus}, y] + 

(1-y) HPL [{plus, minus}, y] +HpL[{Qf pluS/ minus}? y] + 

1 1 \ 

HPL [ {minus, plus, minus}, y] — — HPL [{plus, plus, minus}, y] j 

HPLInt [x / (1 -x A 2) A 3 HPL [{plus, minus}, x] , x] 

x(-7 + 5x 2 ) x(-5 + 3x 2 ) HPL [ {minus} , x] 

32 (-1 + x 2 ) 2 16 (-1 + x 2 ) 2 

7 (31) 

-rT- HPL [ {plus} , x] + \-^pr + —5- HPL [ {plus, minus} , x] 

64 I 32 4 (-1 + x 2 ) 2 J 

• MZV[m] is the Multiple Zeta Value (see Section 2.5) corresponding to the 
index vector m. Their value in terms of mathematical constants are tabu- 
latecCD up to weight 8 and systematically replaced. For higher weights, the 
cases covered by Appendix B are also replaced. 

Table [MZV[ {n} ] , {n, 2, 8}] 

71 -, Zeta[3], — , Zeta[5], -g— , Zeta[7], 



L 5 ' 90 945 9450 

MZV[{-2, -2}] 

tt4 



480 



MZV[{3, 1, 3, 1, 3, 1}] 

tt12 



43589145600 



HPLI [m] is the HPL of weight vector m evaluated at argument i. Their val- 
ues in terms of mathematical constants are tabulated up to weight 8 and 
systematically replaced. For higher weights, the cases covered by Appendix 
A are also replaced. 



HPLI [{plus, minus, plus}] 

n ^ 1 21 Zeta [3] 
-Catalan it + — - — — 
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HPLI[{3, -1}] 

( (-121 + 480 i) 7T 4 + 1440 i 7T 3 Log [2] - 360 tt 2 Log [2] 2 + 360 Log [2] 4 - 
60 i PolyGamma[3, — ] + 8640 PolyLog [ 4, —\ + 3780 (i/r + Log[4]) Zeta[3] 

» The function $HPL0ptions gives a list of the options of the package and 
their current values. 



these tables are those of the FORM package harmpol[68]. 
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$HPLOptions 

{ $HPLAnalyticContinuationSign -» 1, $HPLAutoConvertToKnownFunctions -> False, 
$HPLAutoConvertToSimplerArgument -> False, $HPLAutoLogExtract -» False, 
$HPLAutoProductExpand -> False, $HPLAutoReduceToMinimalSet -> False} 

• The variable $HPLFunctions contains a list of the functions provided by 
the package. 

$HPLFunctions 

{ HPL, HPLlml2pm, HPLAnalyticContinuation, HPLArgTransf orm, HPLAtoM, 
HPLConvertToKnownFunctions, HPLConvertToSimplerArgument, HPLI, HPLInt, 
HPLLogExtract, HPLMtoA, HPLpm21ml, HPLProductExpand, HPLReduceToMinimalSet, MZV} 



3.2 Functions modified 

• We define the derivatives of HPLs as described in Section 2.1. 

{3 X HPL[{2, 10, 4}, x], e x HPL[{-l, 4, -6}, x] , d x HPL[{l, 3, 7}, x] } 

{ HPL[{1, 10, 4}, x] HPL[{4, -6}, x] HPL[{3, 7}, x] } 

I X ' 1+X ' 1 - X ' 

{d x HPL [ {plus, minus, 0}, x] , d x HPL [ {minus, minus, plus}, x] } 

r 2 HPL [{minus, 0}, x] 2 x HPL [ {minus, plus}, x] 

l 1^X2 ' fT^2 I 

The integration showing up in the recursive definition of the HPLs is also 
implemented. It is however recommended to use the function HPLInt. 

r* HPL[{2, 4}, x] r* HPL[{2, 4}, x] r* HPL[{2, 4}, x] 

{ dx, dx, dx} 

l J 1-x J x J 1+x 

{HPL[{1, 2, 4}, t], HPL[{3, 4}, t], HPL[{-1, 2, 4}, t] } 

. r fc 2xHPL[{0, minus}, x] f fc 2 HPL [{plus, 0, minus}, x] , 

U„ 57* J d ^ 

{HPL [ {minus, 0, minus}, t], HPL[{plus, plus, 0, minus}, t] } 

• The function Series is able to expand HPLs around x = and x — 1. 

Series [HPL [{2, 3, 4, 5}, x] , {x, 0, 7}] 

x 4 | 193 x 5 | 1026311 x 6 241822151 x 7 8 
6912 + 1382400 + 8957952000 + 2633637888000 + LxJ 
Series [HPL [{1, 2, 3}, x] , {x, 1, 1}] 

[- \ 4 5 3 3 qq + y HPL[{0}, 0] Zeta[3] + Zeta[3] 2 - -| HPL[{0}, 0] Zeta[5] 
f--^- +HPL[{0}, 0] -HPL[{0}, 0] Zeta[3] + 



y 7T 2 HPL[{0}, 0] Zeta[3] - y HPL [ {0}, 0] Zeta[5] j (1 - x) + 0[1 - x] 2 
Series [HPL [ {plus, minus}, x] , {x, 0, 17}] 

2 x 3 3 x 5 11 x 7 25 x 9 137 x 11 49 x 13 121 x 15 761 x 17 1E 

3 + 5 + 21 + 54 + 330 + 130 + 350 + 2380 + ° LxJ 
Series [HPL [ {plus, minus}, x] , {x, 1, 2}] 

( n 2 3 Log[2] 2 T rol T rl , 1 



+ Log [2 ] Log [ 1 - x] + — Log [ 1 



6 2 ^ l j 2 



x 



2 



^I + yLog[l-x]) (l-x) + (4 + ^I + yLog[l-x ] ) (1 - x) 2 + [ 1 - x] 
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3.3 Working with the options 



The package HPL has some options to control its behavior. They set the pre- 
ferred form in which expressions are displayed. The options can be overridden 
locally by the functions described above. The effects of the options are de- 
scribed in the following. 

• $HPLAutoConvertToKnownFunctions: If set to True, HPLs will be con- 
verted to more common functions (logs, polylogarithms, Nielsen polylogs) 
if possible, using the identities of Appendix C. This might be counterpro- 
ductive when the properties of the HPLs are more explicit in the HPL form 
than other representations. 

$HPLAutoConvertToKnownFunctions = False; 
HPL [{plus}, x] 

HPL[{3, 1, 1}, t] +HPL[{-1, 1}, x] 

HPL [ {plus} , x] 

HPL[{-1, 1}, x] +HPL[{3, 1, 1}, t] 

$HPLAutoConvertToKnownFunctions = True; 
HPL [{plus}, x] 

HPL[{3, 1, 1}, t] +HPL[{-1, 1}, x] 

Log 



1 - x 1 

~T2 + L ° g[ 2 2] Log[2] Log[l + x] +Pol Y Log[2, ^^-] + PolyLog[2, 3, t] 

Furthermore, if this options is set to True while using the analytic continu- 
ation described above, the result may be wrong, as Mathematica does not 
have different conventions for the analytic continuation. Default is False. 
• $HPLAutoProductExpand: If True the products of HPLs are automatically 
converted into a sum of HPL of weight equal to the sum of the weights of 
the two factors, as described in section 2.2. Default is False. 

$HPLAutoProductExpand = False; 
HPL [{plus}, t] HPL [{minus, 0}, t] 

HPL [{plus}, t] HPL [{minus, 0}, t] 
$HPLAutoProductExpand = True; 
HPL [{plus}, t] HPL [{minus, 0}, t] 

HPL[{minus, 0, plus}, t] + HPL [ {minus, plus, 0}, t] + HPL [ {plus, minus, 0}, t] 

Setting the option $HPLAutoProductExpand True affects only explicit prod- 
ucts and does not act on factorized products like the function HPLProductExpand. 



HPL[{-1, 1}, t] (HPL[{1}, t] +t) 

(t + HPL[{1}, t] ) HPL[{-1, 1}, t] 
Expand [ % ] 

tHPL[{-l, 1}, t] +2HPL[{-1, 1, 1}, t] +HPL[{1, -1, 1}, t] 

• $HPLAutoLogExtract: If True the logarithmic divergences log(l — x) and 
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log(rr) are automatically extracted from the HPLs.The default setting is 
False. 



^HPijAutoijogExtract = False; 
HPL[{plus, 0}, t] 
HPL[{1, -2}, t] 




HPL [ {plus, 0} , t] 




HPL[{1, -2}, t] 




$HPLAutoLogExtract = True; 
HPL [{plus, 0}, t] 
HPL[{1, -2}, t] 




HPL[{0}, t] HPL[{plus}, t] -HPL[{0, plus}, 


t] 


HPL[{-2}, t] HPL[{1}, t] -HPL[{-2, 1}, t] - 


HPL[{2, -1}, t] 



The extraction of the divergent behavior only makes sense, if one does not re- 
expand the products automatically with the option $HPLAutoProductExpand 
set to True. If the latter option is set to True, the option $HPLAutoLogExtract 
will have no effect. 



$HPLAutoLogExtract = True; 
$HPLAutoProductExpand = True; 
HPL[{1, -2}, t] 

HPL[{1, -2}, t] 

$HPLAutoProductExpand = False; 
HPL[{1, -2}, t] 

HPL[{-2}, t] HPL[{1}, t] -HPL[{-2, 1}, t] -HPL[{2, -1}, t] 



On the other hand, if the option $HPLAutoConvertToKnownFunctions is set 
to True, the out factorized HPLs of weight one will be replaced by logs be- 
fore being re-expanded, as shown by the following example. 



$HPLAutoLogExtract = True; 




$HPIAutoConvertToKnownFunctions = True; 




HPL[{1, 4, 4}, t] 




HPL [{minus, plus, 0}, t] 




-HPL [{2, 3, 4}, t] -HPL[{3, 2, 4}, t] -2HPL[{4, 1, 4}, t] - HPL [ { 4 


, 2, 3}, t] - 


HPL[{4, 3, 2}, t] -2HPL[{4, 4, 1}, t] -HPL[{4, 4}, t] Log[l-t] 




-HPL[{0, minus, plus}, t] - HPL [ {minus, 0, plus}, t] +HPL[{minus, 


plus} , t] Log [t] 



• $HPLAutoReduceToMinimalSet: If set to True, the HPLs will be automat- 
ically reduced to a minimal basis (up to weight 8). This only makes sense 
if one does not expand the obtained products again, or if the factors of 
smaller weight can be replaced by their expression in terms of known func- 
tions. Therefore, for the reduction to be performed, one has to have the 
option 

$HPLAutoProductExpand equal to False or 
$HPLAutoConvertToKnownFunctions equal to True. 
If this is not fulfilled, the option will have no effect. It defaults to False. 
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$HPLAutoReduceToMinimalSet = False; 
HPL[{plus, 0, minus, plus}, x] 
HPL[{1, 2, 1}, x] 

HPL[{plus, 0, minus, plus}, x] 

HPL[{1, 2, 1}, x] 

$HPLAutoReduceToMinimalSet = True; 
HPL[{plus, 0, minus}, x] 
HPL[{1, 2, 1}, x] 

HPL[{plus}, x] HPL[{0, minus}, x] -HPL[{0, minus, plus}, x] -HPL[{0, plus, minus}, x] 
HPL[{1}, x] HPL[{2, 1}, x] -3HPL[{2, 1, 1}, x] 

• $HPLAutoConvertToSimplerArgument: If set to True, HPLs of arguments 

1/x, x/(x — 1) and (1 — x)/{l + x) will be automatically 
replaced by their representation in term of HPLs of argument x along the 
lines of Section 2.4. Its default is False. 

$HPLAutoConvertToSimplerArgument = False; 

HPL[{1, 2}, 1-x] 

HPL [{plus, 0, minus}, 1/x] 

HPL[{1, 2}, 1-x] 

HPL[{plus, 0, minus}, — ] 

$HPLAutoConvertToSimplerArgument = True; 
HPL[{1, 2}, 1-x] 

1 9 

- — 7T 2 HPL [ {0}, x] - HPL [ {2, 0} , x] - 2 Zeta [3] 

• HPLAnalyticContinuationSign: If set to 1 the analytic continuation of 
the HPLs is taken assuming a positive infinitesimal imaginary part for argu- 
ments, if set to —1 a negative one is assumed. $HPLAnalyticContinuationSign 
is only the default setting and can be overridden by specifying the option 
AnalyticContinuationSign in the function HPLAnalyticContinuation, 

as described above. $HPLAnalyticContinuationSign can be a symbol. The 
default setting is +1. 

$HPLAnalyticContinuationSign 

1 

HPL[{0}, -1] 

i 7i 

$HPLAnalyticContinuationSign = -<5; 
HPL[{0}, -1] 

-i n 6 



3.4 Numerical evaluation 

The package HPL 2 provides an arbitrary-precision numerical evaluation in the 
whole complex plane. 
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3.4-1 Real argument 



HPLs of real arguments with finite precision are automatically evaluated to 
the precision of the argument. 

HPL[{2, 1, 4, 5}, 0.5] 

HPL[{0, plus, minus, 0, minus}, 0.5] 

0.000174907 
0.000809316 

HPL[{2, 1, 4, 5}, — — — ] 

L 2. '30. J 

HPL[{0, plus, minus, 0, minus}, 1 

L 2. '30. J 

0. 00017490 7282457760677740168010076 

0.000809316225590848646 790 73045115 

For arguments with infinite precision, numerical evaluation is only undertaken 
if explicitly requested. 

HPL[{2, 1, 4, 5}, j] 

HPL[{plus, minus, 0, minus}, — ] 

HPL[{2, 1, 4, 5}, y] 
HPL[{plus, minus, 0, minus}, — ] 
N[HPL[{2, 1, 4, 5}, i-], 50] 

n[hpl[{0, plus, minus, 0, minus}, — ] , 50] 

0.00017490728245776067774016801007622686190198808601431 
0.000809316225590848646 790 730451146176 728820505192 76148 



For the numerical values of the MZV, we implemented the procedure described 
in refs. [64,62]. 

HPL[{2, 1, 4, 5}, 1] 
HPL[{0, plus, 0, minus}, 1] 

MZV[{2, 1, 4, 5}] 

121 ^ 4 1 2 T „„roi2 , Log[2] 4 . OTW1 „ T _ r >. 1 



tt z Log[2] z + ^ J + 8 PolyLog 4, — + 7 Log[2] Zeta[3] 



1440 3 " L J 3 1- 2 

HPL[{2, 1, 4, 5}, 1.] 

HPL[{0, plus, 0, minus}, 1.] 

0.0272374 
0.283497 
N[MZV[{4, 3}]] 

0.0851598 

N[MZV[{2, 5, 4, 3}], 40] 

0.0000971639689998491980155438399 78284 78502532 



For real arguments outside the interval [0,1], one can specify the sign of 
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the infinitesimal imaginary part to be used for the analytic continuation 
by appending the option AnalyticContinuationSign as for the function 
HPLAnalyt i cCont inuat i on. 

HPL [{2, 1, 0}, 1.4, AnalyticContinuationSign-* 1] 

HPL [ {minus, plus, 0}, 1.4, AnalyticContinuationSign -> 1] 

-1.99573 - 1.73879 i 
1.08233 - 7.75157 i 

HPL[{2, 1, 0}, 1.4, AnalyticContinuationSign -» -1] 

HPL [ {minus, plus, 0}, 1.4, AnalyticContinuationSign -» <5] //Chop 

-1.99573 + 1.73879 i 
1.08233 - 7.75157 i 5 



3.4-2 Complex argument 

For the sake of the numerical evaluation for complex arguments, we divide the 
complex plane into five different regions as pictured by figure 1. 

region I: \z\ < 0.9 

the series expansion is used 

region II: \z\ > 1.5 

First the argument is brought into region I using the argument transforma- 
tion z — > 1/z described in Section 2.4, then the HPLs are evaluated using 
series expansions. 

region III: 0.9 < \z\ < 1.5 and | arg(z)| < 5tt/12 

First the argument is brought into region I using the argument transfor- 
mation z — > (1 — z)/(l + z) described in Section 2.4, then the HPLs are 
evaluated using series expansions. 

region IV: 0.9 < \z\ < 1.5 and | arg(z)| > 77r/12 

First the argument is brought into region I using the argument transforma- 
tions z — > —z and z — > (1 — z)/{l + z) described in Section 2.4, then the 
HPLs are evaluated using series expansions. 

region V: 0.9 < \z\ < 1.5 and 57r/12 < | arg(z)| < 77r/12 

The evaluation is done using the Holder convolution described in ref. [64]. 

The transformations mentioned above are performed using the "complex" ar- 
gument transformation for the ± weights and the "real" transformations for 
integer weights. The evaluation for the ± weights in region V is done by first 
converting to the integer weights and then using the Holder convolution. 

HPLs with numerical complex arguments with finite precision are evaluated 
automatically. 
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Fig. 1. Complex plane division for the numerical evaluation 



HPL[{1, -2, 3}, 0.5i + 


•7] 


HPL[{plus, minus, 0}, 0. 


51 + 0.7] 


-0.05413 + 0.0438294 i 




0.1566 - 0.508478 i 





HPLs with complex argument with infinite precision are only evaluated if ex- 
plicitly requested. 



HPL[{1, -2, 3}, 1.51 + 1.7] 
HPL[{plus, minus, 0}, 1.51 + 1.7] 

-0.503607 - 0.252553 i 

1. 70143 - 1.07863 i 

N[HPL[{3, -2, 1}, 13 + 71], 50] 

-2. 63 462792 812430377140 92 494101947 7347835545572 74122- 
2.7289184907880926000947141638056212522510887136916 i 

N[HPL[{0, plus, minus, plus}, 13 + 71], 50] 

6.8060017197430759251833294345157620921643906833439- 
6 .49150662339307657806434291734506243502555316220851 



Thanks the possibility of evaluating the HPLs numerically, Plot is able to rep- 
resent HPLs graphically on the real axis 
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Plot[{Re[HPL[{-l, -1, 


-1, 1}, t]], Im[HPL[{- 




i i 11 j_ i n i 
■1, -1, 1}, t] ] }, 


{t, -5, 5}, PlotStyle-* {Hue[l], Dashing[{0 
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and in the complex plane, here are two examples of HPLs plotted on the in- 
terval [—5, 5] and in the unit circle 



- Graphics - 

<< Graphics 1 Graphics ' 

PolarPlot[{Re[HPL[{plus, minus}, 9/10 ® ie ]], Im[HPL[{plus, minus}, 9 / 10 <s 4e ] ] } , 
{6, 0, 2 7T}, PlotStyle -> {Hue[l] , Dashing[ {0 . 02, 0.02}]}] 











-1 -0.5^TV 


V" — 

-\ 


1 — 





- Graphics - 

Plot [ {Re [HPL[ {plus, minus}, ® ie ]], Im[HPL[{plus, minus}, <e 19 ]]}, 
{©, 0, 2 7T}, PlotStyle -> {Hue [1] , Dashing[ {0 . 02, 0.02}]}] 




We checked numerical agreement with ref. [67] at double precision accuracy. 
The values of the MZV have been checked against those of the EZf ace application [69] 
and the GiNaC implementation [62]. Concerning speed, our implementation 
cannot be expected to compete with the GiNaC implementation. 
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4 Conclusion 



In this paper we presented the extension of the existing package HPL to deal 
with complex argument. The new version of the package provides the numeri- 
cal evaluation of HPLs of complex argument over the complete complex plane. 
We also presented new weights more appropriate for complex arguments. They 
will be used for the expansion of hypergeometrical functions around half inte- 
ger parameters [53]. 
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A HPLs at i 



In these tables, we collect the values for the HPLs at i up to weight 3. 



A.l Weight 1 



H(l;i) 



in log(2) in log(4) 



(A.l) 



H{+-i) = -,H(-;i) = -log(2), 



(A.2) 



H(0;i) 



2 



(A.3) 
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A. 2 Weight 2 



H(+,i) = -, ff(-,i)=log(2), (A.4) 

H(+, -i) = i (2C - tt log(2)) , if (-, +, i) = i (-2C + l -ir log(2)) , (A.5) 
if(0,+,i) = 2iC, if(0,-i) = ~, (A.6) 

ff (+, 0, i) = -2iC - ^ , ff (- 0, = " ^ M2) • (A.7) 



A. 3 Weight 3 



H(+,+,+;i 
H(+,+,-;i 

H(+,-,+;i 
H(+,-,-;i 



'48 ' 



= l(7r 2 log(16)-21C(3)) , 



703) 
16 4 



-Ctt + 

7m 3 
" 96 



21C(3) 



2*Clog(2) + ^7rlog 2 (2)-2Li 3 Q-i 



r I 1 » 

+2L l3 (- + - 

7Z7T 3 



(A.8) 
(A.9) 
(A.10) 
(A.H) 

(A.12) 



Y77T' 5 1 i / 1 7 \ 

ff (+, 0; <) = -Cvr + — + -vr 2 log(2) + -«r log 2 (2) + 4Li 3 (- - -) 

(A.13) 



-4Li, 



1 i 

2 + 2 
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H(+,0,+;i) 
H(+,0,-,i) 

H(+,0,0;i) 
H(-,+,+;i) 

H(-,+,-;i) 
H(-,+,0;i) 
H(-,-,+;i) 

H(-,-,-;i) 
H(--0;i) 

H(-,o,+;i) 

H{-A-,i) 
#(-, 0,0;i) 

H(0,+,+;i) 

H(0,+,-;i) 



1 1 77T 3 1 / 

- + 1iC log(2) - log 2 (2) - 2Li 3 ( 
= Cvr , 

= — + 2iClog(2) - -m log 2 (2) + 4Li 3 (- - 

77T 3 1 1 

= Ctt - ^- - Vlog(2) - ^7rlog 2 (2) - 2Li 3 



+2Ll3 b + 2 

:-^ + ^log 2 (2)-2Li 3 C-- 1 - 
96 8 S V ; V2 2 

/ 1 i 
+2L, 3 (- + - 

1 



= -glOg 3 (2), 



^(l2^1og 2 (2)-7r 2 log(4) + 3C(3)) , 
— - r lo g ^(2)-2L„(---) 

/I V \ 



+2Ll3 b + 2 



:^(vr 2 log(2)-3C(3)) , 
:i(7r 3 -6 ? 7r 2 log(2) + 9<(3)) , 



-Ctt + 



7C(3) 



: ~ 32~* ( 7rS + 64C bg(2) + 4?r log2(2) ) ~~ 2Li3 
+2Li 3 



tf(0,+,0;i) = 



1 i 

2 + 2 

.3 



27T 

-Ctt - — 
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H(0, -, (vr 3 + 2vr log 2 (2)) + 4Li 3 (\- % - 



H(0,--;i 
H(0,-,0;i 
H(0,0,+;i 

H(0,0,-,i 
H(0,0,0;i 



8 

-4Li 3 - + - 

3 V2 2 

03) 
16 ' 

in 3 , 3C(3) 
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3C(3) 



27T 

'48 



16 

3 



(A.29) 
(A.30) 
(A.31) 
(A.32) 
(A.33) 
(A.34) 



A. 4 Arbitrary weight 



H( n +;i) = 
H(T-;i) = 
H( n 0;i) = 



{i7i) n 

2 n n\ ' 
(-2 log 2) n 



ni 



2 n n\ ' 

H(»0, +; i) = i2 2n+1 (t(n + 1, 1) - C (n + 



(A.35) 
(A.36) 
(A.37) 
(A.38) 



The HPLs at i with only weights — and can be converted to HPLs at x — — 1 
using the transformation x — > a; 2 . For weight vectors si, s n we have 



if(si,...,s n ;i)=.ff(si,...X;-l) s n = -, s» = 0, -s- = < 



1, Si = - 
0, s 4 = , 
(A.39) 



so long as the last element of the weight vector is not a 0. If it is so, one has 
to extract the 0s from the weight vector using the product identities. 



B Multiple Zeta Values 



We list here some of the identities for MZVs found in [64,65] 
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C(2,1) = C(3), (B.l) 

m-2 

2C(m,l)=mC(m+l)-^2C(m-k)C{k + l), 2 < m e Z , (B.3) 

fc=i 

C(2,"l) = C(n + 2), (B.4) 

C(3, B l) = C(n + 2,l) 

n 4- 2 1 n 

= ^-C(n + 3) - - £ C(* + l)C(n + 2-k), (B.5) 

c( 2) -(2^TT)I UJ ' (R6) 

C(4) -(4nT^)ll2j ' (R7) 

27T 4 ™ 

Cr{3,l})=4-Cr4) = ^- T ^ j , (B.10) 



C(2, 3}) =4" n E(- 1 ) fe C({4n{(4A: + 1)C(4* + 2) 

fc=0 

A; 

-4^C(4i-l)C(4^-4i + 3)}, (B.ll) 



B.l Colored MZVs 



Here we list some identities for MZVs with negative arguments found in [65], 
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2n (_l)n(n+l)/2 

C( n (-2)) = , (B.12) 

sv v ;; 2n + l ! 2 n K 1 



An 

7T 



^"(-- [ )) — ^— ^ ((1 • v'^) +(l-V2) ) 



(4n + 2)! 



_6n 



Cr(-6)) = T^^yyr] (l + 2 3 ^(-l)«(^)/ 2 



I fV3\ 6n+3 + /W3\ 6n+3 _ ! 



(B.13) 



X <I^J + (^J ' (R14) 



C(-i,n) = (-ir 1 fl^, 

n! 

C(-l,-l,"l) = -Li n+2 (l/2), (B.15) 

crc-i)) = (-ir En- 1 ! f " Lifc( [' 1)fc) )" • ( B - 16 ) 

where the sum in the last equation is over all non negative integers satisfying 
Sfc>0 kjk — n 



B.2 Minimal set 



MZV's can be expressed as linear combination of mathematical constants 
(7r, C{ n )-> Li n (l/2)) and, for high weight, a minimal set of other constants. In 
the tables we translated [1], the following choice of constants in the minimal 
set has been used. 



^ = 5(H-i},oo) = E4 L EV L 

il=l ? 1 i 2 =l * 2 

= C(-5, -1) + C(6) ^ 0.98744142640329971377 , (B.17) 

°° r_ 1V 1 h i *2 i 

S7a = 5({-5,i,i},oo) = E L 4 L Ef Ef 

i 1= l H i 2 =l ^2 i3= i «3 

= C(-5, 1,1) + C(-6, 1) + C(-5, 2) + C(-7) (B.18) 
~ -0.95296007575629860341 , 

S75 = 5({5,-l,-l},oo) = EiE^E l ; 



i 1= l 2 1 i 2 =l ^2 i3= i «3 

= C(7) + C(5, 2) + C(-6, -1) + C(5, -1, -1) 

~ 1.02912126296432453422 , (B.19) 
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oo y h y 
s 8a = S({5,3},oo) = E 75 E ~ 



, l= i i i i 2 =i *i 



= C(8) + C(5, 3) ~ 1.0417850291827918834 , 
S8fe = S({-7, -1}, oo) = E £ i^L 

i 1= l *1 ij,=l «2 

= C(8) + C(-7, -1) ^ 0.99644774839783766600 , (B.20) 
s 8c = S({-5,-l,-l,-l},oo) 



»4 



U=l ? 1 ia=l 2 2 is=1 ^3 is=1 «4 

= C(8) + C(-7, -1) + C(-5, -3) + C(6, 2) + C(-5, -1, 2) 

+C(-5, 2, -1) + C(6, -1, -1) + C(-5, -1, -1, -1) (B.21) 
~ 0.98396667382173367094 , 

S8d = 5({-s, -1, 1, i}, 00) = e ^ E 4^ E E - 

il=l 1 i 2 =l ^2 i3=1 «3 j 3= i «4 

= C(8) + C(-5, -3) + C(6, 2) + C(7, 1) + C(-5, -2, 1) + C(-5, -1, 2) 
+C(6, 1, 1) + C(-5, -1, 1, 1) ~ 0.99996261346268344768 . (B.22) 



C Table of representation through more common functions 



In these tables, we collected the identities relating HPLs to more common 
functions that can be found in [70,32]. These identities are used by the func- 
tion HPLConvertToKnownFunctions (see section 3.1)and applied systemati- 
cally when the option $HPLAutoConvertToKnownFunctions is set True, see 
section 3.3. 



C.l Weight 2 



H ({1,1}; x)= l -\og{l-xf 



H({1, -1}; x) = Li 2 ((l - x)/2) - log(2) log(l - x) - Li 2 , 

H ({-1, 1}; x) = Li 2 ((l + x)/2) - log(2) log(l + x) - Li 2 Q , 

/J ({-1, 0}; x) = log(l + x) \og{x) + Li 2 (-x) , 
H{{-2}-x) = -U 2 {-x) . 



(C.l) 
(C.2) 

(C.3) 

(C.4) 
(C.5) 
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C.2 Weight 3 



H ({1, 2}; x) = -2S 1>2 (x) - log (1 - x) Li 2 (x) , 
H({l,l,l};x) 



--log(l-xY , 



H({2,l};x) = S lt2 (x) , 
#({2,-l};x) = Li 3 ( " 



1 + x 



)-Li 3 ( r 



(C.6) 
(C.7) 
(C.8) 



x 



+ x 



Li 3 



1 + x 



-Li 3 (x) + log(l + x)Li 2 (-) + log(l + x)Li 2 (x) 



+ Ilog(2)log(l + x) 2 + Li 3 Q , 
H({-2,l};x) = -S 1 , 2 (x)+U 3 



-2x 



x 



(C.9) 
Li 3 (-x) 



H({-2,-l};x) 
H({l,-l,-l};x) 



+Li 3 Q + Li 3 (x) + log(l - x)Li 2 (-x) 



+ log(l - x)Li 2 i^-j - log(l - x)Li 2 (x) 

+ l - log(2) log(l - x) 2 - 1/6 log(l - x) 3 

■■Si j2 (-x) , 

-i.og(i^).og(l + ^-Li3(i) 



log(l + x)Li 2 



x 



+ Li, 



x 



(CIO) 
(C.11) 

(C.12) 



H ({-1, 1, 1}; x) = 1 log (±±^) log(l - x) 2 + Li 3 (I) 

+ log(l - x)Li 2 (1^) - Li 3 (i^) , (C.13) 

H ({-1, -2}; x) = if ({-2}; x) if ({-1}; x) - 2if ({-2, -1}; x) ,(C14) 
H ({-1, 2}; x) = if ({-1}; x) H ({2}; x) - H ({-2, 1}; x) 

-#({2,-l};x) , (C.15) 
7f({-l,l,-l};x) = if({-l};x)if({l,-l};x)-2if({l,-l,-l};x) . 

(C.16) 



36 



C.3 Weight 4 



H ({l,3};x) = -hi 2 (x) 2 -log(l - x)U 3 (x) , (C.17) 

H({1,1, 2}; x) = ± log (1 - xf Li 2 (x) + 2 log (1 - x) S h2 (x) 

+351,3 (x) , (C18) 

//({l,2,l};x) = -log(l-x)5 li2 (x)-3 ( 5 1)3 (x) , (C.19) 

if ({2, 2}; x) = -2S 2 , 2 (x) + hi 2 (xf , (C.20) 

H({3,l};x) = S 2 , 2 (x) , (C.21) 

if ({2, 1,1}; a;) = S 1)3 (a) , (C.22) 

({1, 1,1,1};*) = ^ log (1-s) 4 • (C.23) 
C.^ Arbitrary weight 

H ({n};x) = U n (x) , n>0, (C.24) 

if ({n}; a;) = -Li_ n (-a;), n<0, (C.25) 

if (f 1}; x) = (-ir log(1 ~ x)W , (C.26) 

n! 

g({ w (-i)};x) = log(1 w j ha;)W , (C.27) 

//({nr 1 !};^^^) • (C28) 
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